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Abstract 

We provide a detailed discussion of the replica approach to thermodynamics 
of a single classical particle placed in a random Gaussian N (>> 1)— dimensional 
potential inside a spherical box of a finite radius L = RyN. Earlier solutions of R = 
oo version of this model were based on applying the Gaussian Variational Ansatz 
(GVA) to the replicated partition function, and revealed a possibility of glassy 
phases at low temperatures. For a general R, we show how to utilize instead the 
underlying rotational symmetry and to arrive to a compact expression for the free 
energy in the limit iV — > oo directly, without any need for intermediate variational 
approximations. This method reveals a striking similarity with the much-studied 
spherical model of spin glasses. Depending on the competition between the radius R 
and the curvature of the parabolic confining potential \i > 0, as well as on the three 
types of disorder - short-ranged, long-ranged, and logarithmic - the phase diagram 
of the system in the (//, T) plane undergoes considerable modifications. In the limit 
of infinite confinement radius our analysis confirms all previous results obtained by 
GVA. The paper has also a considerable pedagogical component by providing an 
extended presentation of technical details which are not always easy to find in the 
existing literature. 

1 Introduction 

In this paper we perform a detailed study of thermodynamics of a single classical particle 
confined to a spherical box filled in with an energy landscape described by a random Gaus- 
sian function TC of iV real variables x = (xi, ...,xn)- Although the problem is meaningful 
for any N, we eventually will be mainly concerned with the limit of large N ^> 1 where 
we will be able to develop a systematic method of analysis. In fact, it is well known that 
such simple yet non-trivial models play a role of a laboratory for developing the methods 
allowing one to deal with problems of statistical mechanics where an interplay between 
thermal fluctuations and those due to quenched disorder is essential. The paradigmatic 
example of systems of this sort are spin glasses pQ, but similar effects are frequently op- 
erational for polymers' behaviour in random environment, for phase separating interfaces 
in random field models or for elastic manifolds pinned by random impurities. In general, 
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presence of quenched disorder leaves no choice but to employ the so-called replica trick, 
which is a heuristic way of extracting the averaged free energy of the system from mo- 
ments of the partition function. The book [2] may serve as a modern introduction to this 
problematic. 

Both conceptual and technical difficulties of dealing with statistical mechanics of dis- 
ordered systems stem from the fact that many features of their low temperature dynamics 
and thermodynamics are dominated by presence of a huge number of metastable states 
(both minima and saddle points of various types ) in the energy functional in configura- 
tion space. At finite temperature those features may generate a complicated free energy 
landscape, and that structure is responsible both for unusual equilibrium properties (an 
"ergodicity breaking", see [TJ [2]) as well as for a complicated long-time dynamical be- 
haviour. The latter manifests itself, in particular, through slow relaxation and aging 
effects [3]. At the level of static properties the broken ergodicity is reflected in an intri- 
cate pattern of spontaneous replica symmetry breaking discovered originally by Parisi [I] 
in the framework of the Sherrington-Kirkpatrick [5] model of spin glass with infinite-range 
interactions. Attempts to relate that picture to the properties of stationary points of the 
corresponding free energy landscape started long ago from the pioneering TAP paper [6] 
and is still a rather active field of research, see e.g. [TJ [8] and references therein for a 
general account, and [HI EDI EH] for a treatment of the model of the present type. 

In this broad context, the model of a single particle in a random potential has played 
an important role in physics of disordered system. It has enjoyed quite a long history of 
research starting from early works by Mezard and Parisi [12], and Engel [13] on static 
properties of such a system, followed by Franz and Mezard [H] and Cugliandolo and Le 
Doussal [3] papers on the corresponding dynamics. To have bona fide thermodynamics one 
has to ensure that the partition function Z of the model is well-defined for any realization 
of the random potential. This is usually achieved by introducing a sort of a confining 
potential V con (x) which prohibits escape of the particle to infinity. The standard choice is 
to use a parabolic potential uniform in all directions V con (x) = |x 2 . The curvature \x > 
then plays together with the temperature T a role of the main control parameter of the 
system, and one of standard goals of the theory is to investigate the phase diagram in the 
(//, T) plane. Technically, the problem amounts to calculating the ensemble average of the 
equilibrium free energy F = —T In Z . The replica trick allows one to represent integer 
moments (Z n ) of the partition function in a form of some multivariable non-Gaussian 
integrals, and one then faces the usual problem of finding ways of evaluating those integrals 
for need of performing the replica limit n — > 0. The crucial step allowing one to achieve 
further progress was suggested by Mezard and Parisi [12] and is widely known as the 
Gaussian Variational Ansatz (GVA). Loosely speaking, it amounts to replacing the non- 
Gaussian integrands with trial Gaussian ones and employing the Feynman-Bogoliubov 
variational procedure to find best possible Gaussian approximation to the true free energy 
in the replica limit, see [121 [33] for a detailed discussion. Proceeding in this way, it turned 
out to be possible to employ again the Parisi scheme of spontaneous replica symmetry 
breaking at the level of the trial free energy, and to reveal the glassy nature of the 
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lower-temperature phase of the model. Since its induction GVA method became one of 
the most popular tools of dealing with quenched disorder. It allowed not only to get 
useful insights into static properties of glassy systems of quite a diverse nature (see e.g. 
[T5J [THl [TTl [THJ, [TH] ) , but was eventually adopted to study the corresponding nonequilibrium 
dynamics as well [3], [201 121] ■ Such a flexibility of GVA is based on the fact that it 
can be formally applied in any dimension N. The accumulated experience shows that 
the physical quantities calculated in such a way usually agree well with the results of 
numerical simulations on qualitative, and sometimes even quantitative level. As a result, 
GVA presently serves as a standard reference point for comparison with any new approach 
to the problem, as e.g. with a recently developed functional renormalization group method 

[22]- 

One should nevertheless recall that the exploitation of GVA lacks formal mathe- 
matical justification for finite dimensions. To this end, Mezard and Parisi provided an 
argument |12j. see also [23], that being a variant of the Hartree-Fock approach the method 
should actually yield exact results in the limit N — > oo. Another intrinsic feature of GVA 
is that the whole procedure is essentially based on a possibility to evaluate explicitly 
some intermediate Gaussian integrals. For the models defined in a restricted geometry 
presence of a geometric confinement may make applications of GVA less convenient. This 
is precisely the case for the class of models to be treated in the present paper: a particle 
confined to an impenetrable spherical box of some finite radius L = RyN, R < oo filled 
with a Gaussian random potential. To make contact with previous works on the prob- 
lem we retain also the parabolic confinement term, so that our model is characterised by 
both R and /i as control parameters. Our main technical observation is that the model 
is exactly solvable for any value of those parameters in the limit N — > oo without any 
need of introducing GVA. Rather, our method is based on observing that the replicated 
partition function possesses a high degree of invariance in the replica space: an arbitrary 
simultaneous O(N) rotation of all n replica vectors x a . An efficient method of dealing 
with such integrals was developed recently in the work of one of the authors [2H [2S] 
within the framework of the theory of random matrices. This method allows us to arrive 
to the effective free energy functional in the replica space in the most economic, and 
we believe elegant way. The subsequent analysis follows the standard route of using the 
Parisi Ansatz for the spontaneous replica symmetry breaking pattern, and the nature 
of the low-temperature phase is known to depend very much on the decay properties of 
the covariance function of the disorder potential[T21 13]. Actually, we propose a simple 
mathematical criterion which allows one to formally discriminate between the long-ranged 
and short-ranged disorder, and also suggests to single out potentials with logarithmically 
growing correlations as a separate intermediate class. For all types of the disorder we 
provide a detailed derivation of the free energy functional, the stability analysis, and a 
thorough description of the most important features of the resulting phase diagram. The 
latter undergoes considerable modifications reflecting a competition between the confine- 
ment provided by the radius R and the curvature of the parabolic potential \i. For the 
limiting case of infinite confinement radius our results faithfully reproduce all those fol- 
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lowing from earlier applications of GVA. This is just another explicit verification of the 
expected exactness of GVA in the limit N — > oo. 

One more point we find appropriate to mention here is that the Parisi scenario of 
replica symmetry breaking for spin glass models is changing presently its status from a 
powerful heuristic method of theoretical physics to an essentially rigorous mathematical 
procedure, well-controlled in the case of models of infinite range. This important change 
is mainly due to recent seminal results by Talagrand [26, 27J, based on earlier works by 
Guerra [28], see also interesting works by Aizenman, Sims and Starr [29]. In particular, 
Talagrand was able to demonstrate that the equilibrium free energy emerging naturally in 
the Parisi scheme of replica symmetry breaking is indeed the correct thermodynamic limit 
of the free energy of the paradigmatic Sherrington-Kirkpatrick model and, more recently, 
for the so-called spherical model of spin glasses, originally studied by one of the present 
authors in [30]. It is natural to expect that similar justification should be possible also 
for other types of models with glassy thermodynamics, the model of a particle in random 
potential being the most natural candidate. However, the solution of the problem in the 
framework of GVA does not seem to be a promising starting point for such a verification. 
On the other hand, the expression for the free energy emerging in our approach is actually 
very close in its form to those emerging in the spherical model of spin glasses (see discussion 
in the end of the next section). The revealed similarities of our problem to the spherical 
spin glass model give strong evidence in favour of applicability of Talagrand's method for 
the model under consideration. We leave a detailed investigation of this issue for a future 
work. 

Finally, we hope that our presentation has also a certain pedagogical value by providing 
extended description of a few technical details known to experts, but which are not always 
easy to find in the existing literature. 

2 Definition of the model and its formal treatment 
by replica method. 

As was discussed in the introduction, we consider a classical particle confined to an im- 
penetrable spherical box of some finite size L. To ensure the non-trivial behaviour in 
the limit N —>■ oo, one has to scale the radius of the sphere with N, and we denote the 
corresponding domain as {Djy : x 2 < N R 2 } [§. 

As usual, the main object of interest for us is to calculate the ensemble average of the 
free energy 

F = -T\nZ, Z=[ exp-/?ft({x})dx, (1) 

Jd n 

where (3 = 1/T stands for the inverse temperature. The average of the logarithm over the 
disorder (denoted in the present paper by angular brackets) is performed with the help of 

2 For some applications, as e.g. [10] and [11], it is useful to keep in mind that the volume Vn(R) of 
such a sphere in the large N limit behaves asymptotically as Vn ~ where L e — y/2-Ke R. 



4 



the standard replica trick, i.e. the formal identity 

i „ n 

(InZ) = lim - In (Z n ), Z n = / e^^-i w{Xa} TT dx a . (2) 
^on Jd n a=l 

The standard choice of the energy function for this problem is 

H{x} = jE4 + V(x ll ...,x N ) (3) 
z k=i 

with fi > 0. A random Gaussian-distributed potential V(x) is characterized by zero mean 
and the variance specified by the pair correlation function which we choose in the form 
ensuring stationarity: 

(y(x 1 )F(x 2 )) = iV/(^(x 1 -x 2 ) 2 ) . (4) 

The previous analysis (T2J [13] revealed that one should essentially distinguish between 
two rather different situations. The first describes the case of a short-range correlated 
disorder corresponding to functions f(x) vanishing at infinity, with f(x) = e~ x being a 
typical representative of that class. In the second case the correlations are long ranged, 
and at large distances the potential grows in such a way that: 

([V(x 1 )-V(x 2 )] 2 )cx(x 1 -x 2 ) 2 ^, (5) 

the exponent 7 to be chosen in the range < 7 < 1. The particular case 7 = 1/2 
corresponds to a potential V(x) being the standard Brownian motion. Although formally 
V(x) in the latter case can not satisfy the property Eq.(J3]) as its variance (V 2 (x)) is 
obviously position-dependent, one can easily satisfy oneself that such a difference is com- 
pletely immaterial for the free energy calculations. As a result, we always assume for the 
long-ranged disorder the validity of Eq. (jlj) with the choice: 

f(x) = f(0)-g 2 x\ /(0)>0, 0< T <1. (6) 

In what follows we will be able to formulate a certain criterium relating the nature of 
the low-temperature glassy phase of the model with the shape of the correlation functions 
f(x), see Eq. (l59l) and the discussion around it. According to that criterium, in a broad 
class of random potentials with short-range correlations the glassy phase will be described 
by the so-called one-step replica symmetry breaking (1RSB), whereas all long-ranged 
potentials Eq^S]) are characterized by the full replica symmetry breaking (FRSB). The 
criterium also suggests naturally to single-out as a special case the logarithmically growing 
correlations, that is 

f(x) = -g 2 \n(x + a 2 ), a 2 < 1. (7) 

We shall see that such a choice leads to the phase diagram which combines some features 
typical for the short-ranged behaviour and others for the long-ranged types of disorder. In 
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this sense the logarithmically growing correlations should be considered as a marginal case 
intermediate between the two broad classes described above. It is interesting to mention 
that the choice of a glassy model with logarithmically growing correlations Eq.([7|) is not 
purely academic, as such objects actually emerge e.g. in a context of statistics of the wave 
function in disordered two-dimensional systems |31j . 

Actually, a few initial steps in evaluation of the free energy are precisely the same 
for all types of disorder, provided the latter is of the Gaussian nature. Performing the 
averaging over the disorder in Eq.flS]), we in the standard way arrive at the following 
expression: 

a 2 r n 

(Z n ) = e WW / exp -{3H n {x a } [J dx a , (8) 

JD N a =l 

where 

^n{x a } = |E X a-^E/(^(Xa-X 6 ) 2 ) . (9) 

a=l a<b 

So far all our manupulations were exact. To achieve further progress one has to suggest 
an efficient way of working with the resulting multidimensional integral. In the stan- 
dard model with infinite box radius R = oo Mezard and Parisi[T2] suggested to deal 
with apparently non-Gaussian character of the integrand by replacing the exact repli- 
cated Hamiltonian if n {x a } with a trial Hamiltonian if^{x a } chosen to be Gaussian with 
respect to all variables x a , and then to apply a kind of variational principle to find the 
best possible Gaussian approximation. 

Here we point out the possibility of a different route, valid for any value of the param- 
eter R and requiring at this step no approximation. It is based on observing that the inte- 
grand in Eq.flSJ) in fact possesses a high degree of invariance: it depends on N— component 
vectors x a only via n{n + l)/2 scalar products q a b = x a Xb, a <b, and is therefore invari- 
ant with respect to an arbitrary simultaneous O(N) rotation of all vectors x a . Moreover, 
our choice of the integration domain respects this invariance. 

An efficient method of dealing with such integrals is based on the possibility of rewrit- 
ing the integral 



J Nn = ... Z E (xi,...,x n )dxi...t2x n (10) 

J\xi\<L J\x„\<L 

whose integrand X r (x 1 , ...,x n ) possesses such type of invariance in an alternative form as 

J N ,n = C N , n [ J Q (Q) detQ( N - n ^ 2 dQ, (11) 

provided N > n + 1. Here in Eq. fllll) the original integrand Xj.({x a })) is expressed as 
a function Xq(Q) of n x n real symmetric positive semidefinite matrix Q > whose 
entries are precisely those g a b, introduced above. The integration domain D ® is simply 
= {Q — O5 laa < NR 2 , a = 1, . . . n}, the volume element is dQ = n a <bdq a b and the 
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proportionality constant is given explicitly by 



7T 2 



N- 



'N,n 



nn—1 p / N-k \ 
k=0 1 V 2 ^ 



'12) 



Although some use of similar formulae in statistical mechanics can be traced back to 
[22] , the full potential of the transformation from Eq. ffTUl) to Eq. ffTTT) seems to be revealed 
in the work of one of the authors [2H [25] where it was rediscovered in the context of 
random matrix theory using Ingham-Siegel matrix integrals. Since then the formula 
found applications in physics of disordered systems, see e.g [33] for an elegant derivation 
and further use. 

Applying this procedure to our case and using the subsequent rescaling Q — > NQ 
yields the following exact expression for the averaged replicated partition function: 



(Z n ) = C Nn N Nn ' 2 e^ Nn ^ [ (detQ)- {n+1)/2 e -/w*»«) dQ 



where 



*n{Q) = £ E Qaa - 7^ In (det Q)-(3j2f 



a=l 



2(3 



a<b 



-{q aa + qbb) - qab 



(13) 



(14) 



and TV is assumed to satisfy the constraint TV > n. The final integration domain Dq is 
already iV— independent: Dq = {Q > 0, q aa < R 2 , a — 1, . . . n}. 

The form of the integrand in Eq. (1131) is precisely one required for the possibility of 
evaluating the replicated partition function in the limit iV — » oo by the Laplace ("saddle- 
point") methoco Taking into account the expressions Eqs.([T|), ([2D, and f[T2l) the free 
energy of our model is then given by 



lim Lf) = -^ln(27re) - ^/(0) + hm -$ n (Q) 



(15) 



where the entries of the matrix Q are chosen to satisfy the conditions: = for 

a < b. This yields, in general, the system of n(n + l)/2 equations: 



1 



Q~ 



-(qaa + qbb) - qab 



o. 



1,2, 



n 



(16) 



3 One should note that as long as one is interested only in finding the leading exponential factors 
in TV — * oo limit, one can in principle follow a different route. Namely, impose n(n + l)/2 constraints 
Qab = ^x a Xb through the integral Fourier representations involving n(n + l)/2 auxiliary fields X a b, and 
take the saddle-point both in A— and Q— variables. The same asymptotic result Eq. ([T5)l then follows 
after simple manipulations. We however believe that the use of the elegant mathematical procedure based 
on Eq. pTj) is more conceptually clear, and appealing aesthetically. We hope that it may also provide a 
useful basis for calculating 1/TV corrections, and perhaps for a rigorous mathematical treatment of the 
problem. 
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and 



! Q- 1 L + f3f 







ab 



^{Qaa + Qbb) - Qab 



0, a^b (17) 



where f'(x) stands for the derivative df Jdx. 

One should also ensure that the solutions to these equations respects the constraint 
Qaa < R 2 for all a = 1, . . . , n imposed by presence of the boundaries of the integration 
domain Dq. We will shortly see that in the replica limit n — > that condition will be 
violated for some regions of parameters T, fi. If this happens, Eqs. (I16j) should be simply 
replaced by equalities q aa = R 2 , and the solution to the remaining set Eq.( fl7|) should be 
sought with those constraints imposed. One then notices that for R = 1 the resulting 
expression for the free energy formally coincides , up to a trivial constant term, with that 
obtained for the so-called "spherical" mean field model of spin glasses. Various features of 
the latter model attracted a lot of research interest in recent years, see e.g. [30j EH ESI [27] . 

The correspondence between the free energies of the two models which is so apparent in 
our approach deserves a short comment. In one of its recent incarnations the Hamiltonian 
H(cr) of the spherical model was defined as [2"7] 

n ( a ) = E N (p-i)/2 E 9h,..,ipVh---<ri P (18) 

P>1 h,—,ip 

in terms of iV— component vectors a spanning the sphere of radius yN. Here gi lr .,,i p 
denote independent Gaussian variables with mean zero and unit variance. The original 
version [30] of the same model can be shown to yield precisely the same free energy in 
the thermodynamic limit [27]. Immediate consequence of the definition Eq. (fl8l) is that the 
correlation function (H(cri)7Y(cr 2 )) g depends on the vectors <J\ and <7 2 only via the scalar 
product (o"i<T2). This is precisely the property sufficient to ensure global O(N) invariance 
of the replicated partition function of the model, and our method of deriving the mean 
free energy goes through without any modifications, and may well be the shortest possible. 
The spherical constraint J2i<n °f = N simply translates into the condition R — 1. 



3 Analysis of the phase diagram of the model within 
the replica symmetric ansatz. 

Our procedure of investigating the equations Eqs. (1 16 [ 171) in the replica limit n — >• will 
follow the standard pattern suggested by developments in spin glass theory. We first 
seek for the so-called "replica symmetric" solution, and then investigate its stability in 
the (//, T) plane. When the replica symmetric solution is found inadequate, it will be 
replaced by the hierarchical ("Parisi", or " ultrametric" ) ansatz for the matrix elements 
q a b, with various levels of replica symmetry breaking. Since the full analysis contains a 
lot of features to be explained in detail, the reader may wish to cast a regular look at the 
resulting phase diagrams figl, fig. 2 and fig. 3 in the process of reading. 
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The Replica Symmetric Ansatz amounts to searching for a solution to Eqs. fll6|ll7l) 
within subspace of matrices Q such that q aa = q d , for any a = 1, . . .n, and q a< b = go, 
subject to the constraints < g < q d < R 2 . Inverting such a matrix Q yields again the 
matrix of the same structure, with the diagonal entries all given by 

* + »(»- 2) 
(Qd- qo){qd + qo{n-l)) 

and off-diagonal entries given by 

Po = —, rj—^. - f 7TT (20) 

{qd-qo){qd + qo{n-l)) 

Note, that 

Pd~Po = (21) 

Qd - qo 

The system Eqs.f ll6|fT7l) is reduced in this way to two equations, which we write directly 
in the replica limit n — > as 

p-Tp d + ^f(q d -q ) = (22) 
-Tp + ^f'(q d -q ) = (23) 



This system of equations is easy to solve employing the relation (BID , anc ^ to obtain 

In order this solution to be sensible one first of all has to require f'(x) < 0, which we 
always assume to hold in our model. In addition, the solution Eq. fT24l) can hold only as 
long as q d < R 2 , which in view of the above expressions amounts to the condition 



T 1 „, T 



R 2 - - + -J' - > . (25) 

When the above inequality is violated, we should rather use q d = R 2 , and find go from 
the "spherical model" type equation 

(/ " + 7^ f (R 2 ~ lo) = (26) 



(i? 2 -g ) 2 T 2 

following immediately from Eqs. ()23|l20p . 

Let us briefly discuss general properties of the boundary line T 6 (/x) separating the 
"/i— dominated" regime from the " R-dominated" one in (/i, T) plane, for a fixed value of 
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R 2 -- = -- 2 r(-) • (27) 



the confining radius R. For a given value of \i the value of TJ, is obtained by solving the 
equation 

t jl = _Lf (It 

When analysing this equation we shall assume in addition to the condition f'{x) < 
two more conditions: the uniform concavity condition f"(x) > as well as the condition 
f"'(x) < 0. Sensibility of that choice will be justified by analysis of patterns of spontaneous 
replica symmetry breaking in subsequent sections, see discussions around Eq. (1591) . We 
also assume that f'(x) — > for large x, as is indeed the case for all types of disorder in 
the original model. 

For the short-range disorder the values /'(0) < oo and /"(0) < oo. It is convenient to 
define for subsequent use two quantities 



A*e = 4V-/'(°)> R 



™ (28) 



Then a simple graphical analysis of Eq. 1 )271) shows that for fi > \i c that equation has a 
single solution T&(/j) tending asymptotically to Tb = fiR 2 for /x ^> \i c . In contrast, for 
fx < fi c the number of solutions essentially depends on the value of the confining radius 
R. Defining R cr as in Eq. (l28l) . we find that for R < R cr the condition /j < /j c implies that 
no such solution Tb exists at all. This fact corresponds to the picture of monotonically 
increasing curve T&(/x) starting from the point (/i c , 0) in the ip,T) plane. In the opposite 
case R > R cr there are two solutions T&i < Tb2 in the whole interval /iq < fi < /i c , with 
Tb2 — T b \ — ► as /i — > /i - The value of /i and the corresponding temperature value 
T = Tfci = T 62 can be found as 



Vo = \/f"(r ), T = /i to, (29) 
with t being the solution of the equation 

R 2 = h(r ), h{r)=r-§p- (30) 

Note, that dh/dr = / / ( r )/ / "( r )/[/"( r o))] 2 > according to our assumptions. Thus, the 
right-hand side of Eq. (!30|) is a monotonously increasing function, and thus the equation 
has a (unique) solution r > as long as R 2 > h(0) = R 2 r . In contrast, for R < R cr the 
equation Eq. (l30~l) has no solutions. 

Finally, for fi < /i the equation ( [271) has no more solutions. 

Relegating similar analysis of the long-ranged as well as the logarithmic correlations 
to the end of the section, we now discuss the last important ingredient of the procedure: 
the issue of the stability of the emerging solution against fluctuations in the replica space. 
Indeed, the very essence of the saddle-point method calls for a check of the replica sym- 
metric solution being locally stable, in the sense of corresponding to the true extremum 
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(in the replica limit n — > 0, to a maximum) of the functional This should be judged 

by analysing the eigenvalues of the matrix G a b tC d = gJ^§~ ' a ^ ^ ms solution. Such analysis 
along the lines of the classical De-Almeida and Thouless (AT) paper [36] is presented in 
the Appendix B. The main outcome is that the /i-dominated replica symmetric solution 
Eq. (|24p is locally stable as long as 

f-r(l)>o, (si) 



whereas the "spherical model"-type solution satisfying Eq. ff26|) is stable provided 

' " T^f" (R 2 - ?o) > 0. (32) 



(i? 2 - g ) 2 T 2 

Combining the latter with Eq. (l26|) for q one easily finds that the domain of the replica- 
symmetric stability in the R-dominated regime is given by 



T > T yJf»(T ) (33) 



where To is precisely the solution of Eq. (j30|) . 

We will see shortly that both conditions Eq. fl3TT) and Eq fl33|) could be violated for low 
enough values of /i and T below the so-called de-Almeida-Thouless line Tat(/j)- Actually, 
in the R— dominated region of the phase diagram such a line is parallel to the //—axis, 
as the solution r of Eq. ( 130]) determining the right-hand side of ( 1331) depends only on R 
and hence is \l— independent. In particular, for the limiting case \i = the expression 
for the whole de-Almeida-Thouless line in (R,T) plane is just Tat(R) = T o\f f"( T o) (see 
e.g fig.4(c) for a particular choice of the short-ranged potential). 

Let us see how these features are incorporated into the analysis of the boundary T&(/x) 
between the /^-dominated and the R— dominated regions performed by us above. In the 
case of a short-range disorder with /"(0) < oo studied above we can easily see that the de- 
Almeida Thouless line Tat(^) given by the equality sign in Eq. (l3~Tj) must end up for zero 



temperature at the point hat — yf"(ty- As a simple consequence, for small confining 
radius R < R cr , with R cr again given by Eq. (l28|) . the whole instability region falls outside 
the domain of validity of /^-dominated solution. To find whether the replica-symmetric 
solution is stable one therefore has to use instead the equation given by the equality sign 
in Eq.( l32|) . As the latter equation actually does not have at all roots for R < R cr , we 
conclude that for the present model the replica symmetric solution is always stable for 
such values of the confining radius R. In short, the small box size R < R cr precludes 
possibility of the replica symmetry breaking, hence the glassy behaviour of the models 
with short-ranged correlations, see fig.3c@ 



4 Note that in recent works [TUJ, [TT] precisely the same lengthscale R cr appeared in the analysis of the 
"geometric complexity" associated with the zero-temperature limit of the present model. Namely, the 
samples with R < R cr were found to be unable to support the existence of exponentially many saddle 
points in their energy landscape. Such an existence seems to be a necessary feature of a phase with glassy 
behaviour [IT]. 
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The situation changes crucially for R > R cr . First of all, now the lines T),(/i) and 
the AT line Tat{(J>) can indeed intersect each other in (T,fj) plane. Curiously enough, 
the point of intersection turns out to coincide with the point (fi ,T ), see Eq. (f2Ul) . which 
emerged in the previous analysis as the "leftmost" point on the curve T{,(/z), such that 
for < /i the equation (T27j) has no solutions. For fi > /x the branch T&i drops therefore 
below the AT line and hence has no meaning. At the same time for another branch 
Tb2 > Tat and hence that solution survives. We shall see in subsequent sections that 
our analysis of the phase with broken replica symmetry will provide us with a modified 
expression for the boundary line T 61 (/i) extending from the point (/x ,T ) down to zero 
temperature. 

For \x < fiQ we are in the domain of validity of the R-dominated solution, and the cor- 
responding AT temperature should be given by \i— independent value Tat from Eq. (!33|) . 
The latter expression now makes full sense as the corresponding solution To indeed exists. 
Finally, in view of Eq. (1291) it is evident that the \x— independent value of Tat is simply T 
everywhere in R-dominated phase, so that the R-dominated and \x— dominated AT lines 
indeed meet each other at the same point (//o, Tq) of the phase diagram, as was natural 
to expect. 

So far our analysis assumed the case of a short-range disorder. For the case of long- 
ranged disorder, Eq.(j5j), the overall structure of the line is in many respects similar, 
but has some peculiarities. In particular, there will be no analogue of the critical value 
of the confinement radius R = R cr in this case. Below we proceed in presenting a brief 
analysis of the situation for any value of the exponent 7 G (0, 1). Using f'(x) = — r yg 2 x 1 ~ l , 
the equation Eq. (l25l) for the boundary line 7&(//) between the \i— dominated and spherical 
model like regions in the (/i, T) diagram can be written as 



F(u) = vJ^ - R 2 u + 7^ = 0, with u = [ — ) (34) 




We see that the function F(u) has its single minimum for u > at u = u m i n = 
(fz 2 -R 2 ) 7 ; an d T(0) > 0, F{u — ► 00) > 0. Hence, the equation F(u) = has two 
positive solutions «i )2 as long as F{u min ) < 0, and no solutions if F{u min ) > 0. Then a 
simple calculation shows that two branches T&i < T^ exist as long as 

^7 1/2 (2-7) 1 - 7/2 

^ > ^ = r,9 7-, 7^ ■ ( 35 > 

R 2 -~i(l--i) — 

At /j, = Hq those two branches merge: T bl = T b2 = T which implies u 2 = U\ = u min . The 
corresponding characteristic temperature can be found as 

rp _ ^7 1/2 (2-7)~ 7/2 ^ ^ 

(I-7) 2 

Finally, we should take the presence of the de Almeida-Thouless conditions Eqs.f l3"Tf3"31 
into consideration. Recall that for the case of a long-ranged disorder f"(0) = 00 in 
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contrast to the finite value typical for the short-range case. As a consequence, the AT 
line now extends in the (/x, T) plane to arbitrary large values of /x and is explicitly given 
by the equation 

7 

T AT = ^ (37) 

[^2 7 ( 1 _ 7 )] 7 - 2 

It is now easy to verify that the boundary line T b (ii) which is given for /i > fio by 
the two branches T b i and T b2 intersects AT line precisely at the point (/i ,T ), with 
uat = (^at/a*) 1-7 - Moreover, one can check that F{uat) < for /i > /x . This simply 
implies that T b i < Tat < T b2 for any /x > /x - Hence, the lower branch T b i should 
be discarded as falling outside the domain of validity of the replica symmetric solution. 
The analysis of the phase with broken replica symmetry presented in the subsequent 
sections of this paper reveals that the line T b i(n) shall be rather replaced by a vertical 
line n — hq extending from the point (/xo,To) down to zero temperature. Finally, for 
all values /x < fio the de- Almeida - Thouless line Tat{h) is given by the /x— independent 
temperature Tat = T (R). In particular, this expression just provides the de-Almeida- 
Thouless line in (R, T) plane for the limiting case /i = 0. 

At last, we provide the analysis of the replica symmetric solution for the logarithmic 
correlations, Eq.([7|). In this case both AT line and the boundaries T b i,T b2 can be easily 
found explicitly in the whole /x— dominated regime. They are given by 



R 2 -a 2 ± J(R 2 + a 2 ) 2 - V//x 2 



(3* 



where the first formula holds for /x < fx at — g/a 2 , and in the second formula the upper 
sign corresponds to T b2 , and the lower to T&i. The two branches of T&(/x) meet at fi = 
Ho = 2g/(R 2 + ° 2 ) so that there is no T b solution for // < /i . The de-Almeida-Thouless 
line meets the boundary T b (fi) a ^ precisely /io- In order to make such intersection happen 
one has to ensure that /io < hat, which is possible only if the confinement radius R 
exceeds the critical value R cr = a. As long as R > R cr for \i > /j one can see that 
T fel (/i) < Tat(h) < T b 2 (h), so that only the upper branch T b2 makes actually sense. A 
subsequent analysis of the phase with broken replica symmetry will again reveal that the 
line T b \(fi) should be replaced by the vertical line /j = /i everywhere in the glassy phase. 
Finally, for \i < Ho the de-Almeida-Thouless temperature is given by the \i— independent 
value Tat — T — g-$E^?- Again, for the limiting case [L — this expression provides the 
de-Almeida-Thouless line in the whole (R, T) plane. In particular, for R —>■ R cr = a we 
have To —>■ 0, showing that there is no place for broken replica symmetry in the logarithmic 
case as long as R < R cr . 

To summarize our findings, we present in fig. 1, fig.2 and fig. 3 the resulting phase 
diagrams in (/x, T) plane for the particular choice (i) f(x) = /(0) — g 2 x 1 ^ 2 of long-range 
disorder, (ii) for the logarithmic correlations, and finally (iii) f(x) = e~ x corresponding 
to the short range disorder. In all these cases the AT line can be found explicitly. For 
the /i— dominated regime in case (i) Tat = J'/J^/* anc ^ ^ n case ^at = — 2/x In \i- 
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(b) 

T 




R 



(c) 

Figure 1: The phase diagrams for the long-ranged potential with 7 = 1/2 and g = 2. 
The case (a) corresponds to R = a/3, the case (b) to R — oo, and the case (c) to the 
choice /i = 0. Dotted line in the case (a) represents the wrong branch of the boundary 
Tb between the R— dominated and /x— dominated phases with broken replica symmetry, 
and is replaced by the vertical full line. The notation for phases are as follows: I stands 
for fi— dominated replica-symmetric (RS) phase; II for //— dominated glassy phase with 
broken RS; III for R— dominated RS phase, and IV for R— dominated glassy phase with 
broken RS; 



14 



The corresponding coordinates of the intersection point (ij, ,T ) for the two models: (i) 

= T o = 0^ and (iii) /x = e"^ 2 " 1 )/ 2 , T = (R 2 - 1) e"^ 2 " 1 )/ 2 . Note, that the 

branch Tn of the solution to (l2Tj) has no particular meaning as it drops as a whole below 
the AT line. In the presented phase diagrams we also included the replacement of that 
line by the correct boundaries T fe (/i) everywhere in the glassy region T < Tat- 

We also presented for completeness the typical phase diagrams for systems with long- 
range and logarithmic correlations in two limiting cases: in the (//, T) plane for R = oo 
and in (R, T) plane for \x = 0. The de-Almeida-Thouless temeperature in the latter 
case is equal to Tat — T Q (R) = 9 ^/a 2 for the long-ranged potential with 7 = 1/2. The 
corresponding diagram for the case of short-range correlations will be presented in the 
end of section Sec. (14.21) . 

We finish this section with writing down explicit expressions for the equilibrium free 
energy in the replica-symmetric solution. They are given correspondingly by 

foo = -~ In (2nT/n) - ^ [/(0) - /(T///)] (39) 

for the fi— dominated region; in particular, Foo^o = /'(0)/2/x. As for the R— dominated 
RS region, the free energy is given by 

^00 = \r\h - \) - T - In (2ntT) - ± [/(0) - f(tT)] , (40) 
with t being the solution of the equation ( see Eq. ( 127)1) ) 

^ - j + f\tT) = (41) 



In particular, t\ T=0 = Rj J — /'(0) and Foo|t=o = \^R 2 — Ry—f'(0). Along the line T ft (/i) 
given by Eq.( |27l) we obviously have t — and the two expressions for the free energy 
indeed coincide. In particular, the two expressions for Foo\t=q are indeed equal for /i = fi c , 
see Eq.(j2SD- 



4 Free energy functional within the Parisi scheme of 
replica symmetry breaking. 

In the region below the AT lines one must discard the unstable replica symmetric solution 
in favour of one with the broken symmetry. In the present section we derive the expression 
for the free energy of our model in the phase with broken replica symmetry. We will 
follow a particular heuristic scheme of the replica symmetry breaking proposed originally 
by Parisi [1] in the framework of the SK model and employed for the case of the spherical 
spin glass in [30] . 
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In our dealing with the free energy functional Eq. (fl4l we follow closely the method 
suggested in the paper by Crisanti and Sommers [30]. To make the present paper self- 
contained we choose to describe the procedure in Appendix A and to provide a few 
technical details skipped in [30] • The analysis is based on explicit calculation of the 
eigenvalues of the Parisi matrix Q. 

We are actually interested in the replica limit n —>■ 0. According to the Parisi prescrip- 
tion explained in detail in the Appendix A, in such a limit the system is characterized 
by a non-decreasing function of the variable q denoted as x(q). Such a function depends 
non-trivially on its argument in the interval q < q < qk, with q > and qk < qa- Outside 
that interval the function stays constant: 

x(q < go) = 0, and x(q > q k ) = 1. (42) 

In general, the function x(q) also depends on the increasing sequence of k positive 
parameters m 8 satisfying the following inequalities 

< rrii < m 2 < . . . < m k < m k +i = 1 • (43) 

As shown in the Appendix A, in the replica limit the following identity must hold for 
any differentiable function g(x): 



lim -Tr[g(Q)) 

n-+0 n 



9 {Qd -qk)+ 9 



'ik 



x(q) dq I dq 



(44) 



In particular, for the first two terms entering the replica free energy functional Eq. (1141) 
application of the rule Eq. (144 p gives 



lim — 

n->0 n 



TrQ - —Tr In (Q) 



2P 



'ik 



i 



tz dq ■ (45) 



2/T" v ™ w ' 2(3 Jo f q qd x(q)dq 

The last term in Eq. (fl4l is also easily dealt with in the Parisi scheme (see Appendix A), 
where it can be written as 

k 



lim f 

ajtb 



Maa + qbb) ~ qab 



i™E(%i ~ m i)f(qd 



qi 



1=0 



'Id 







f(Qd ~ q)x'(q) dq, 
(46) 

by using explicitly the derivative of the generalized function Eq. (11091) . Using integration 
by parts and taking into account the properties Eq. (]42l we finally arrive at the required 
free energy functional for the phase with broken replica symmetry 



ix T 
2 Qd ~2 



ln(2ne(q d - q k )) 



o qd ~ qk + In" x (q) dq 



- dq 



+ ^ ( f(Qd 



q k ) - /(0) 



<ik 



f{qd-q)x{q) dq. 



(47) 



This functional should be now extremised with respect to the non-negative non- decreasing 
continuous function x(q). In the \x— dominated regime it also should be extremised with 
respect to q&, whereas for the R— dominated situation the latter variable is fixed to be 
q d = R 2 . 
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4.1 Solution within the scheme of fully broken replica symme- 
try. 

The actual analysis essentially depends on whether we consider the number of steps k 
in the Parisi scheme of the replica symmetry breaking (see Appendix A) to be finite, or 
allow it to be infinite. In the latter situation one conventionally speaks about a case of full 
replica symmetry breaking (FRSB), and in the rest of the present section we restrict our 
consideration to this specific case. The number of steps k is assumed to tend to infinity 
in such a way, that the sequence of the parameters mi is replaced with a continuous 
variable u G [0,1], such that mi + \ — m\ = Ami — * du. The sequence of parameters 
qi satisfying Eq. (j!05j) is simultaneously transformed to a non-decreasing differentiable 
function q(u) changing between its minimal value q{u = 0) = go an d its maximal value 
q(u = 1) = qk- The function x(q) defined in Eg. (11 09 p is then assumed to be transformed 
to a non- decreasing function continuous in the interval go < g < qk, and at least once 
differentiable there. Outside that interval it satisfies Eq. fj42|) . As is easy to see, the 
identity a;[g(u)] = u holds as long as g < q{u) < qk- 

Variation of the free energy functional Eq. fl4"7"l) with respect to such a function x(q) as 
well as with respect to qd yields 



SF. = \ 



Qk 

5x(q)S(q)dq + Sq d P 

10 



0, (48) 



where 

S(q) =T f ^ - 2 + h\q d - q) (49) 

[qd ~ Qk + Xf x(q) dq\ 

and 

T fik dq 1 1 fik 

P = V + T / 7 qk +-f {q d -qk) + 7^ / x(q)f (q d -q)dq. 

qd-qk Jo [q d - q k + }q ik x(q) dq\ 2 T T J qo 

(50) 

Stationarity therefore always amounts to the condition 

S(q) =0, VgG [q ,q k ] , (51) 

and for the \i— dominated regime we must add also the equation P = 0. As S(q) = 
implies ^S(q) = we can differentiate Eq. fj49|) once, and get after a simple algebra the 
equation 

rih _ T 
q<i~qk+ / x(q) dq = Vg G [g , q k ) ■ (52) 

h s/Fiqd - q) 

Next differentiation immediately yields the explicit formula for the function x(q) 

rp f"'(q d — g) 

x ® = ~2 [f"{q d -q)]W Vg G qk] ' (53) 
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What remains to be determined are the values for parameters qo, qj., q^, the last one only 
for the fx— dominated case. To this end, we substitute the value q = qu into the relation 
Eg. (152]) and obtain a simple equation for the difference q^ — q^. 

T 

qd-Qk= 7= • (54) 

y/f"(qd - ?*) 

Next, we use the condition S(qo) = 0. Exploitation of Eqs. (l49ll52l) yields after a simple 
algebra one more relation: 

_ fig* - Qo) r ,,x 
f' iid - go) 

Finally, we can use already obtained relations to simplify drastically a somewhat cum- 
bersome equation P = 0, see Eq. (l50|) . which holds only in the //-dominated regime. In 
particular, the explicit form Eq. (1531) for the function x(q) allows one to write: 

^ £ x{q)f{q d ~ q)dq = -\ £ = ^to " ^ ~ [T (?d " %)f 2 

Now taking into account the condition S{q^) = and the relation Eq. (1541) allows one to 
see that the stationarity condition P = amounts simply to 

^ = f'ild ~ qo) • (56) 
which in turn allows to rewrite Eq. fl55|) in the fi— dominated regime as 



qo = — \f{qd - qo) • (57) 

The system of three equations Eq. fl54l . Eq. fl56l) and Eq. (j5"7|) for the \x— dominated regime 
or of the two equations Eq. fl54l . Eq. fl55|) for the R— dominated regime allows one to 
determine the parameters qo,qk (and, if necessary, q d ) as long as the function f(x) is 
specified. The corresponding equilibrium free energy can be written, after some algebraic 
manipulations exploiting Eq. ()52ll53ll54|l55l) . in the form valid for both regimes: 

Foe = 2 qd ~2 ln t 27re ( gd ~ qk ^ + 2T ~ Qk ^ ~ ^ ~ ^ d ~~ qk ^'^ qd ~ qk ^ 

-qo^f"{qd-qo) - r yjf"(q d ~ q) dq, (58) 
J go 

where for the R— dominated phase we assume that qa = R 2 . 
This seemingly completes our solution. 

Consistency of the Parisi FRSB scheme requires, however, that the emerging function 
x(q) given by Eq. (l53l must be a real- valued non-negative non- decreasing function of its 
argument. Non- negativity and reality of this function is indeed ensured by earlier imposed 
conditions f"(x) > and f"'(x) < 0. Non-negativity of the derivative of the right-hand 
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side in Eg. (1551) yields as the necessary consistency condition a new ineguality A(q d —q) < 0, 
where the function A(x) is expressed in terms of f(x) as 

A{x) = \[f"{x)] 2 -f\x)r{x)- (59) 

Checking our three main choices for the function f(x) - the long-ranged, the short ranged 
and the logarithmic - we observe that the above ineguality is indeed strictly satisfied for 
any long-ranged potential Eg. ([6]) where 

A(x) = ^V(7 - 1) 2 (7 - 2)x 2 ^ 3 ) < 0, V 7 € (0, 1) . 

The ineguality is only marginally satisfied in the logarithmic case Eg. (17]) where it is easy 
to see that A(x) = 0. And it is strictly violated in a typical short-range potential, as 
e.g. for f(x) = e~ x when A(x) = \e~ 2x > 0. We therefore conclude that the full replica 
symmetry breaking can not occur in short-range models, neither in the fi— dominated, 
nor in R— dominated regimes, and should therefore be replaced by a different scheme. 
We will discuss the necessary modifications in the next section. As to the long-ranged 
potentials, the scheme is fully legitimate and solutions to the above eguations can be 
easily found for any value of 7 G (0, 1). In the \l— dominated regime they are given by 
following expressions: 

\vh ( fi 2 \^ / T 2 \^ . , 

( g ^)— (i_ 7 )^2 \g 2 i{i-i)j 7(1-7)/ 



and 



. , T 2-7 1 . . 

X \Q.) = 77~T7, 71 7711 7T, < q < Qk ■ (61) 



2 [ 5 r2 T ( 1 _ 7 )] 5 (q d -q) 
A few more useful relations follow from those above 



g d = (2-7)g , x(q ) = [g 2 -f(l - -f)}~ , and x(q k ) = 1 - - . (62) 

l ^7-2 l 



According to the general procedure such a solution makes sense as long as q^ > q , and it 
is easy to check that the condition can be rewritten as T < Tat, where Tat is precisely 
the de-Almeida-Thouless temperature for this model given by Eq. (j37p . 

A curious feature of the above solution is that the values of q , x{q^) and q d turn out to 
be temperature-independent everywhere in the phase with broken replica symmetry. As q d 
is nothing else but the thermodynamic expectation value of the collective displacement in 
the original model: q d = ^x 2 ^> , it is conventional to speak in this case about "freezing" 
of the system below the AT line. On the other hand, the same feature ensures that in 
this "frozen" or glassy phase the boundary line T&(/x) in the (fi,T) plane separating the 
jj,— dominated regime from the R— dominated one must be vertical. Indeed, it is given by 
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the condition q d = R 2 which together with Eq. (l60l) simply amounts to /i = /x , with the 
value of fiQ given by Eq. (l35l) . Everywhere in the R— dominated glassy phase we have 



. . T 2-7 1 . . 

x{ - q > = ~oTTT, 7^2 ^' 90 - q - qk ' ( 63 ) 

2 [0 2 7(l-7)]a (R 2 -q)i 



with 



2 



R 2 ^ T~ 



Qo — 7. ? Qk — R — 



2-7 ' [^ 7(1 _ 7)] - 

These observations complete our investigation of the phase diagram for the case of long- 
ranged correlations, full version of which was presented for the particular case 7 = 1/2 in 
fig. la. 

Now we consider the case of logarithmic correlations, Eq.([7j) assuming R > a. As for 
such a potential A(x) = this implies ^x(q) = for q < q < qk- Indeed, our general 
relations Eqs. (J5S]) . ([51|) . (J5B|) and Eq. (j57|) give in this case an especially simple solution 

T g g a 2 g 2 

x(q) = — , q <q< qk] where q = -; qk = 2 ^7— and q d = 2 a . (64) 

/i fj, l-T/g \i 

for the /i— dominated regime. The consistency condition qk > go is satisfied as long as 
T < Tat, where Tat is as expected the de-Almeida-Thouless temperature. For this case 
Tat is given by Eq. fl38|) . The condition q d = R 2 which defines the boundary T b (/i) with 
the R— dominated glassy regime again amounts to /i = /i = 2g/(R 2 + a 2 ), with the same 
value of /io we found earlier for this case. Finally, everywhere in the R— dominated glassy 
phase occupying the rectangle T <T = g ^ 2 ~° 2 , fi < /i in the (fi,T) plane we have: 

T , # 2 + a 2 2 2 T/g 

x{q) = — , qo<q< q^ where g = = , = R ~a - — — . (65) 

9 2 l-T/g 

Using these expressions one can find the corresponding equilibrium free energy F^. The 
latter takes especially simple form for zero temperature: 

■2^2 for fi> hat = 9/a 2 , 
? oo|t = o = { -fa 2 -(?ln(^) for /i < /i < fi A T , (66) 
§ i? 2 - </ [l + In i^M-ti] for < /i < /i = 2g/{R 2 + a 2 ) . 

Here the upper line corresponds to the \x— dominated RS regime, second line to the 
\i— dominated solution with broken RS symmetry, and third line to R— dominated so- 
lution with broken RS, respectively. 

The corresponding phase diagram was presented in fig. lb. 
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4.2 Short-ranged correlated potentials: solution within the one- 
step replica symmetry breaking. 

In this section we will discuss the solution of our problem pertinent for the potentials 
with the correlation function f(x) corresponding to the uniformly positive values of the 
function A(x) defined in Eq.( l59|) . As we argued above, this situation is typical for the 
short-ranged correlated potentials. 

We have seen that searching for a solution within the FRSB scheme when A(x) > 
leads to a contradiction. The only remaining possibility within the Parisi hierarchical 
ansatz is therefore to assume that the number k of steps in the replica breaking hierarchy 
is finite: 1 < k < oo. We shall see, adopting for our model the line of reasoning suggested 
first by Crisanti and Sommers in |30j, that the condition A(x) > forces us to select 
the value k = 1 as the only possible. This situation is conventionally called the one- 
step replica symmetry breaking (1RSB). In the Appendix B we also sketch the stability 
analysis showing that 1RSB solution is indeed stable versus small fluctuations in the full 
replica space. 

In the case of a general finite k the function x(q) is again a non- negative and non- 
decreasing, but has only a finite number of points of growth in the interval [^o^fe], 
parametrised by the sequences of qi and mj £ [0,1]- The variation of the free energy 
with respect to those parameters is again given formally by the same formula Eq. (148]) but 
with understanding that 

k k 

Sx(q) = ^2(8mi +1 - 6mi)8(q - q{) - XX m '+i ~~ m 0<K<? - Qi)^Qi ■ (67) 

1=0 1=0 

The corresponding part of the variation of the free energy is therefore given by 

1 k rik 1 k 

SF=- J2($ m l+i - $ m l) / S(q)dq- - J2( m l+i ~ m l) S(q^)Sqi , 
1 i=o Jqi 1 1=0 

where the function S(q) is defined in Eq.( l49l) . As Stuq = Srrik+i = 0, the above expression 
can be represented in the form 

26F = ]T 6m t / S(q) dq - £(m m - m,) S( qi )6 qi , (68) 
i=i J i 1 - 1 1=0 

showing that the stationarity conditions amount to the system of equations 

/ S{q) dq = for I = 1, 2, . . . , k - 1, and 3{q{) = for I = 0, 1, . . . , k . (69) 

As the function S(q) is obviously continuous in [qo,qk], the first type of condition en- 
sures that it takes both positive and negative values in each of the intervals [qi-\,qi\. 
According to the second condition S(q) vanishes at both ends of those intervals, therefore 
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this function must have at least one maximum and at least one minimum in each of the 
above intervals. Denoting position of those extrema as q e , we must have j-S(q)\ qe = for 
each point of extremum. Then differentiating Eq. (|49p we get after a simple algebra the 
equation whose solutions should give us all possible positions of extrema g e : 

[i* T 
qd~qk+ x(q)dq = Vg e e go,9* • (70) 

J *> s/FWd - q e ) 

Now, it is convenient to consider both the right- and left-hand side of this relation as some 
functions of the variable q e . Taking the derivative over q e from the left-hand side once 
gives — x(q e ) < 0, showing that the left-hand side is a decreasing function of its argument. 
Next differentiation yields —-^-x(q e ) < 0, which shows that the left hand side is a concave 
decreasing function. Now we treat the right-hand side of Eq. ( 1701 in the same way. First 
differentiation shows that the right-hand side is also a decreasing function of q e due to 
the condition f"'{x). At the same time its second derivative turns out to be equal to 

T 1 

- A(q d - q e ) — - > 

2 [f"(qd - q e r /2 

due to the condition A(x) > 0. Thus we see that the right-hand side is a convex decreas- 
ing function of q e . As any convex decreasing function can have at most two points of 
intersection with a concave decreasing one, there is no more than two different solutions 
of the Eq. tJTD]) for q e . Hence, there could be only a single interval [go, gi] in the description 
of the function x(q). Using henceforth the notation mi = m, we see that such a function 
is given simply by 

f 0, q < g 

x(q) = < m, g G [g , gi] . (71) 

{ i, q > qi 

This is precisely the one-step replica symmetry breaking (1RSB) hierarchical ansatz, which 
is thus shown to be the only possibility within the Parisi scheme for systems satisfying 
A(x) > 0. 

In the rest of this section we study the phase diagram resulting from the implemen- 
tation of such a 1RSB scheme in our model. Instead of using the general stationarity 
conditions Eq. (1691) . we prefer to start directly with the variational free energy functional 
depending on the parameters go, gi, qa, and m G [0,1]. In fact, we find it more convenient 
to use the set of variables 

m; g ; Q = q% - g ; and y = q d -qi 



as independent variational parameters. Substituting Eq. (j7T|) into Eq. (j47| and using the 
above notations we arrive at the variational free energy of the form 



F^ bB = ^(q + y + Q)--\n(2ne)--(l--) \ny -— In (y + mQ) 
' (/ " +^l(l-rn)f(y)+mf(y + Q)-f(0)). (72) 



2 y + mQ 2T 
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One should also remember that in the \i— dominated regime the found solution should 
respect the inequality 

qo + y + Q = q d <R 2 . (73) 

In the R— dominated regime the above inequality must be replaced with the equality, and 
the arising constraint reduces the number of independent variational parameters to three. 

Let us start first with the fi— dominated regime. The parameter go enters in the 
functional linearly, and can be immediately excluded in favour of the relation 



T 

y = y m {Q) = — rnQ 
n 



(74) 



This fact allows one to write down the variational free energy as the function of only two 
variables, m and Q. It is also natural to operate with the so-called excess free energy given 
by the difference AF = F^ SB — F^f between the free energy value of the RS solution, 
Eq.( l39l) and the variational 1RSB free energy, Eq.( }T2l . After simple manupulations we 
obtain 
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mQj + mf I - 4 ( I 
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(75) 



Note that for m — 1 the difference AF vanishes, since this choice obviously brings us 
back to the replica-symmetric solution. Differentiating Eq.( J75l) over Q and assuming 
(1 — m) 7^ we obtain the first stationarity condition 



Q 



y m {Q) 



(f[yo(Q)}-f \y m {Q)\) 



(76) 



where y m {Q) is the combination from the right-hand side of Eq. (174l) . and we introduced 
also the notation yo(Q) = y m (Q) + Q- Similarly differentiation of Eq.(J75l) over m yields 
another equation 



T i 

2 ln 



^y m {Q) 



Q {m + f f [y ° m ) " T (/ [y ° m " ; [Vm{Q)]) ■ (7?) 



This set of equations determines the system behaviour in the fi— dominated regime of the 
phase with one-step broken replica symmetry. Finally, to find the transition line 71 (/i) 
to R— dominated regime one needs to check the condition Eq. f[T3"j) . For this one should 
be able to express the value of go i n terms of Q and m. To this end we notice that the 
stationarity condition with respect to Q at the level of the original free energy expression 
Eq. (172|) together with Eq. (!74l) immediately produces the required relation: 



1 



- 2 f'[yo(Q)} 



(7E 
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Before discussing features of the resulting phase diagram in full generality it is worth 
pointing out the existence of a particular case when the system of equations Eqs. (176|77j) 
and Eq. (T74"]) allows for an explicit algebraic solution. This is precisely the case of log- 
arithmically correlated potential, Eq.([7j). Indeed, a direct verification shows that the 
substitution 



Q = 9 - 



m = — 



T 

9 



(79) 



i-W 

solves all those equations. This observation should look less surprising if one comes 



back to the solution Eq. flMl) found earlier as a particular limiting case of the full replica 
symmetry breaking ansatz. One then realizes that the corresponding function x(q) was 
constant in the interval q G [qo,qk], the feature being characteristic rather of one-step 
replica symmetry breaking. This is just another evidence towards marginal nature of 
the potential with logarithmic correlations: the case can be looked at both as a limiting 
special case of FRSB solution, and that of 1RSB solution. 

After this digression, we proceed to discussing typical features of the phase diagram 
in a generic case of one-step replica symmetry breaking. As we have already seen, the 
transition to the phase with broken replica symmetry may occur in two different ways: 
either along the line Q — > 0, or along the line m — > 1. Let us first investigate the behaviour 
of the free energy difference f l75|) for small Q < 1. Expanding that expression to up to 
the first two nonvanishing terms gives 
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Q 2 



2T 



-/i 2 + /"(T//x) 
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-m^ + {\-m)f"{T^) 



(80) 



Maximization with respect to Q shows that Q = for /i 2 > /"(T//i) which is precisely 
the de-Almeida-Thouless condition Eq. (l3Tj) of stability of the replica-symmetric solution. 
Below the AT line the maximum of the excess free energy happens at a nonzero value of 
Q given by: 



;±-m)f"(T» 
This allows one to write AF in terms of m only as 
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m(l — m) 



y + f"(T/v)f 
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(82) 



We see that the excess free energy grows cubically in the glassy phase in the vicinity of 
the AT line. Such a behaviour is typical for continuous glass transitions to a phase with 
broken replica symmetry. 

Requiring maximum of this excess free energy with respect to m one finds the equi- 
librium value of this parameter in the vicinity of the AT line: 



m 




(83) 
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Such value of m may tend to the maximal possible value m = 1 when the temperature 
T and the parameter /i approach along the AT line a point (/x m ,T m ) where they satisfy 
the system of two equations: 

= 0, -/4 + r(— 1=0. (84) 

It is easy to check that these conditions are precisely those ensuring that (/x m , T m ) is the 
point of a maximum of the AT line, i.e. dT J^ = 0, d2 J^ T < 0. Close to the point of 
maximum the AT line is described by 

35 2 





S > 0, (85) 



where t A t = (T m - T AT )/T m < 1, 8 = (fx - fi m )/Hm < 1- 

For the point of maximum to happen within the \x— dominated regime the correspond- 
ing value fi m should obviously satisfy \i m < liq, with /x being the value earlier defined by 
Eqs. ( 129(501) . This condition after a simple algebra reduces to the inequality 

r 2 > Rl = -2^4 " §r\ ^ 

f"{n) f"[r m ) 

for the confinement radius R. Here r* is the (unique) solution of the equation t* = 
— 2 jjjrPr\ = ^(t*). Existence (and uniqueness) of such a solution is ensured for the short- 
ranged potentials satisfying our main condition A(x) > 0. Indeed, in that case the 
function z{r) satisfies z(0) > and dz/dr = 1 — 2v4(t)/[/" / (t)] 2 < 1, so that the point of 
intersection of the graph of the function y — z(t) and the straight line y — r is obviously 
unique. 

As long as R cr < R < R* the position of the maximum on the AT line is irrelevant as 
the AT line earlier meets the line T&(/x) at the point (/j,q,Tq). As we have already seen in 
our analysis of the replica-symmetric solution, the line T&(/i) must have at the meeting 
point the vertical tangent: j-T B (n)\ m = oo. The line T b (fi) then continues to the phase 
with broken replica symmetry. The corresponding equation is given by solving the system 
of equations Eq. flToT) and (1771) together with the condition: 

y (Q)- - 2 f[y»{Q)] = R 2 - (87) 

The resulting transition line T&(/x) evidently ends up at zero temperature T = at some 
point fib of the fi— axis. A typical phase diagram of this sort was presented in fig.lc for the 
particular case of f(x) = e~ x , when the value R* can be found analytically as R* = v^3, 
see Eq. (l86|) . Note that in contrast to the long-ranged case now the line Tb(p) is not 
strictly vertical inside the glassy phase, although actual numerical values of its derivative 
are quite big. Let us also note that in the R— dominated part of the glassy phase (i.e. to 
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the left of the line T b (/z)) the values of the parameters m, Q, q at a given temperature T 
"freeze" to their \i— independent values which those parameters acquired precisely at the 
point of the transition line T b (/i) such that T b = T. 

It is appropriate to recall that the AT stability line Eq. (1311) for a short-range potential 
meets T = line at the point \x = [±at = [f'ify] 1 ^ 2 - Therefore for such systems a 
non-trivial transition to the phase with broken replica symmetry is possible also at zero 
temperature with decreasing /i, in contrast to the case of long-ranged potentials. In 
fact, it is known that some aspects of zero-temperature behaviour may be amenable to 
investigation without resorting to replicas, see [HI ED], so it may be useful to provide below 
a more detailed picture of the zero-temperature transition within the replica approach. In 
performing the zero-temperature limit one first should realize that \im.T^om/T = v < oo, 
so that the relevant parameters governing the system behaviour in that case will be v and 
Q. The variational excess (free) energy is given in terms of these parameters by 

2AF T=0 = nQ + ^hy(l-nvQ)+v [f(Q) - /(0) - Qf'(0)] . (88) 

The corresponding stationarity conditions can be reduced after some algebra to the set 
of two equations: 

7 ^ ? , = / / (Q)-f(0), (89) 
1 — fiv Q 

l - In (1 - n v Q) = -fi Q + v [f(Q) - /(0) - Qf'(Q)} , (90) 

which also could be obtained by performing the limit T — > directly in Eq. (176|) and ( [77]) . 
A straightforward calculation shows that close to the transition point /i = hat = [/"(0)] 1 / 2 
the parameter v tends to the value v = — { -P\ and the corresponding excess energy at 

^ AT 

zero temperature behaves again cubically in the vicinity of transition (cf. Eq. (1821) : 

3/"'(0) \HAT J 

Let us now consider the case of infinite confinement radius R — > oo, which was the 
subject of earlier studies within GVA [I3] . In this situation it is interesting to investigate 
the behaviour of the system at zero temperature and vanishing confining potential \i — > 0. 
Analysis of the system Eq. (189j) then shows that as long as fi — > necessarily Q — * oo in 
such a way that fivQ — > 1, provided for large Q holds f(Q) — > as well as Qf'(Q) — > 0. 
One then finds after some algebra that the parameter v must satisfy in this limit the 
equation 

., r ~i(o)v 2 _ a 

f(0)' 

The valid solution v — 3> 1 is given asymptotically by 

<-jwA-m>) +0{,nl, ^ h (92) 
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This in turn yields for the parameter Q and for the excess free energy the following 
asymptotic expressions: 



/(o) m V f'(o)J 

Recalling that for zero temperature qa = Q — \f'(Q), hence Q — > oo implies qa ~ Q, 
we indeed see that Eq. (!93|) agrees with one obtained in the framework of GVA[13] for the 
mean displacement parameter = ~ (x 2 ). One also can perform similar calculations in 
the opposite limiting case \x = and R — > oo, the result for T = being given by the 
same asymptotic formulae Eq. (1931) with the only replacement fi 1/R 2 . 

Returning to considerations of general finite values of the parameter R, we see that 
the \x— dominated regime described at T = by the equations Eq. fl89l) can exist as long 
as 

\f(Q) < P 2 

where we took into account the zero-temperature limit of Eq. flTB]) . With decreasing \x the 
latter condition is first violated at some value = which can be found numerically by 
solving the system of three equations: Eq.( l89|90l) and ( |94l) . For lower values of //, i.e. for 
ji < fi b , the parameters Q and v will "freeze" at their boundary values v = and 

Q = Q{^b)- 

The shape of the phase diagram described above as typical for R < R*, see Eq.( )86|) . 
experiences essential modifications as long as R exceeds R*. In the latter case the AT 
Temperature line Tat{^) starts to decrease with decreasing \x to the left of the maximum: 
jj, < fi m . We shall see that the correct form of the phase diagram in this situation requires 
considering the transition line to the phase with broken replica symmetry given by the 
condition m — > 1. 

To investigate such a possibility we expand the excess free energy Eq. flToT) in powers 
of p = 1 — m as 

AF = A(Q) p + B(Q) p 2 + C(Q) p 3 + . . . , (95) 

Extremising with respect to p gives the condition A(Q) + 2pB(Q) + 3p 2 C(Q) + . . . = so 
that the equilibrium value of p for small p is approximately given by p m ~ — A(Q)/2B(Q) + 
. . .. The excess free energy in the phase with broken replica symmetry is then given by 
AF w A(Q)p rn + B(Q)p 2 n = -A(Q) 2 /4B(Q) > 0, which implies B(q) < 0. Asp m > 0, we 
conclude A(Q) > 0. It is evident also that the line at which p m — > implies the condition 
A(Q) = 0. Moreover, extremum conditions of Eq. (!95|) with respect to Q implies in the 
limit p m — > another condition A'(Q) = 0. This consideration shows that the line at 
which m — > 1 within the phase with broken replica symmetry is given by the system of 
equations: 

A(Q) = 0, ±A(Q) = 0, (96) 



27 



dAF I 



where explicit expression for A(Q) can be found from the condition A(Q) — -^-| p=0 . 
This gives, up to an unessential factor 



f T -Q 
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f - +Qf - 



(97) 



and therefore 
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t 2 i - m 
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r( T - 



(98) 



Substituting these expressions to equations Eq. fl96|) we see that the resulting conditions 
are indeed equivalent to the pair of equations (l76l) and ( |77l) . provided m — > 1. We will be 
looking for its non- vanishing solution Q = > 0. 

The equation Eq.flSSD for Q -> reads A'{Q) « $ [-/i 2 + f"(T/fi)], hence is 
positive just below the AT instability line Eq. (1311) . We may expect this sign to change 
in the vicinity of the point of maximum (fi m ,T m ) from Eq. (lM]) . as we know that the line 
m = 1 passes through that point. In the vicinity of AT line the parameter Q is small, 
and we can expand 
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One indeed finds from this expression a non-vanishing solution Q* > of the system 
Eq.(l96l) for T < T m , fi < \i m close to the point of maximum to be 



o 2# + ir(5) 



. tun ( %n\ _ 



r 2 



(100) 



Substituting this back to Eq. (!97|) we get an explicit relation between T and /i defining 
the line Ti(/x) along which a solution with broken replica symmetry Q ^ appears first 
with m = 1 close to (/x m , T m ): 
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(101) 



Expanding in powers of T\ = (T m — Ti)/T m <C 1, 5 = (// m — p)/fi m ^ 1 one can find to 
the leading order the equation for this line to be 



6 2 



T 2 (T N 



6/i 



AV, 



5 < 0, 



(102) 



which shows that the transition line T\{p) when approaching the point (fi m , T m ) from the 
left has in that point its maximum. At the point of maximum it therefore smoothly meets 
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the AT line which in the same approximation is described by Eg. (1851) . with only second 
derivative experiencing a jump. 

For lower values of /i away from \i c the temperature decreases down to a point 

where that line meets the boundary line and the system transits to R— dominated 

regime. The point of intersection can be found only numerically. For even lower values of 
\i the transition temperature Ti(p) freezes to its value at the intersection point, see fig.4. 
In the limiting case of infinite confinement radius R = oo the fi— dominated regime covers 
the whole phase diagram, and it makes sense to investigate the asymptotic behaviour of 
the line Ti ([/,) for fi — > 0. Anticipating the behaviour Til^o — > in such a way that 
T\ » fi one can show that the corresponding no n- vanishing solution > of the 
system Eq. (l96l) must be sought close to its maximal possible value = Q™ ax = T/fj,. 
Using /(oo) = f'(oc) = one finds after some algebra the corresponding asymptotic 
expression for the transition line T = Ti(/x): 

f(0) /(0) T 
P= ^exp — — , T->0. (103) 

For completeness, we present in fig.4 also typical phase diagrams for systems with 
short-range correlations in two limiting cases: in (T, fj,) plane for R = oo and in (T, R) 
plane for fi = 0. Recall, that in the latter case the de Almeida-Thouless temperature 
for our choice of the correlation function is explicitly given by Tat — Tq(R) = (R 2 — 
l)e-!(* 2 -D. 
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A Appendix: Parisi matrix, its eigenvalues and eval- 
uation of traces in the replica limit. 

We start with describing the well-known structure of the n x n matrix Q in the Parisi 
parametrisation. At the beginning we set n diagonal entries q aa all to the same value 
Qaa = 0. This value will be maintained at every but last step of the recursion. The 
off-diagonal part of the matrix Q in the Parisi scheme is built recursively as follows. 
At the first step we single out from the n x n matrix Q the chain consisting of n/mi 
blocks of the size m\ < n, each situated on the main diagonal. All off-diagonal entries 
Qa/3, a 7^ P inside those blocks are filled in with the same value q a p — q\ < 0, whereas 
all the remaining n 2 (l — l/m±) entries of the matrix Q are set to the value < < q\. 
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The latter entries remain from now on intact to the end of the procedure, whereas some 
entries inside the diagonal m x x mi blocks will be subject to a further modification. At 
the next step of iteration in each of those diagonal blocks of the size m 1 we single out the 
chain of m 2 /mi smaller blocks of the size m 2 < m 1; each situated on the main diagonal. 
All off-diagonal entries q aj 3, a 7^ (3 inside those sub-blocks are filled in with the same value 
Qap — 12 > Qi, whereas all the remaining entries of the matrix Q hold their old values. 
At the next step only some entries inside diagonal blocks of the size 7712 will be modified., 
etc. Iterating this procedure step by step one obtains after k steps a hierarchically built 
structure characterized by the sequence of integers 

n = m > nil > m 2 > • • • > n^k — m k+i — 1 (104) 

and the values placed in the diagonal blocks of the Q matrix satisfying: 

< g < qi < Q2 < ■ ■ ■ < q k (105) 

Finally, we complete the procedure by filling in the n diagonal entries q aa of the matrix 
Q with one and the same value q aa = qd> qk- 

For the subsequent analysis we need the eigenvalues of the Parisi matrix Q. Those can 
be found easily together with the corresponding eigenvectors built according to a recursive 
procedure which uses the sequence Eq. (11041) . It is convenient to visualize eigenvectors as 
being " strings" of n boxes numbered from 1 to n, with I th component being a content of 
the box number I. 

At the first step i — 1 we choose the eigenvector to have all n boxes filled with the 
same content equal to unity. The corresponding eigenvalue is non-degenerate and equal 
to 



Ai = q d + qk{fn k - 1) + g fc _i(m fe _i - m k ) + . . . + <?i(mi - m 2 ) + qo{m - mi) (106) 

Now, at the subsequent steps i — 2, 3, . . . , k + 2 one builds eigenvectors by the following 
procedure. The string of n boxes of an eigenvector belonging to i th family are subdivided 
into n/m.i_i substrings of the length m^i, and numbered accordingly by the index j = 
1,2, ... , n/rrii-i. All mj_i boxes of the first substring j = 1 are filled invariably with all 
components equal to 1. Next we fill rrii-i boxes in one (and only one) of the remaining 

— 1 substrings with all components equal to —1. In doing so we however impose 

a constraint that the substrings with the indices j given by j = 1 + l^f^ 1 should be 

excluded from the procedure, with / being any integer satisfying 1 < I < — 1. After 

the choice of a particular substring is made, we fill all n — 2mj_i boxes of the remaining 
substrings with identically zero components. It is easy to see that all di = n/mi_i—n/mi_2 
different eigenvectors of i th family built in such a way correspond to one and the same 
d— degenerate eigenvalue 

\ = Qd + qk{fn k - 1) + g fe -i(m fe _i - m k ) + ... + q i -. 1 (m i - 1 - m,) - ^-2(^-1) (107) 
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In this way we find all n possible eigenvalues, the last being equal to 

Afc+2 = Qd-Qk fn k+1 = q d - q k . (108) 

The completeness of the procedure follows from the fact that sum of all the degeneracies 
di is equal to 

n ^\ ( n n \ ( n n s 



1+ 1 + +•••+ = n 

\mi J \m 2 m\J \m k +i m k J 

Note that all the found eigenvalues are positive due to inequalities Eq. (1105|) between 
various which is required by the positive defmiteness of the matrix Q. Note also 
that all eigenvectors built in this way are obviously linearly independent, although the 
eigenvectors belonging to the same family are not orthogonal. The latter fact however 
does not have any bearing for our considerations. 

To facilitate the subsequent treatment it is convenient to introduce the following (gen- 
eralized) function of the variable q: 

k 

x (q) = n + J2(m l+1 -m l )6(q-qi) (109) 

1=0 

where we use the notation 9{z) for the Heaviside step function: 6{z) = 1 for z > and zero 
otherwise. In view of the inequalities Eq. fll04ll05j) the function x(q) is piecewise-constant 
non-increasing, and changes between n and 1 as follows: 

X {<1 < Qo) =m = n, x(q < q < qi) = m 1 , . . . , x(g fc _i < q < q k ) = m k , x(q > q k ) = m k+1 = 1 

(110) 

Comparison of this form with Eq. (jl09j) makes evident the validity of a useful inversion 
formula: 

1 1 k ( 1 1 \ 

^r = - + E (in) 

which will be exploited by us shortly. 

As observed by Crisanti and Sommers[30] one can represent the eigenvalues Eq. (11071) 
of the Parisi matrix in a compact form via the following remarkable identities: 

rid rid rid 

Ai = / x(q)dq = nq + / x(q)dq, \ i+2 = \ x(q) dq, i = 0,l,...,k (112) 

JO Jqo Jqi 

As a consequence, these relations imply for any analytic function g(x) the identity 

-y -y k+2 y , „g d k / y 1 \ / fid \ 

-Tr[g(Q)] = -^g(\ i )d i = -g (nq Q + x(q)dq)+J2\ )g( x(q) dq) . 

n n i=i n \ Jq / l=Q \rri[ + i mi j \J qi j 

(113) 
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Next one observes that taking the derivative of the generalized function from Eq. (Ill II) 
produces 



d 



x(q) 



;ii4) 



This fact allows one to rewrite the sum in Eq. ( 11 13ft in terms of an integral, yielding 

1 



1 1 / fid \ fQk+0 / fid \ d 

-Tr [g{Q)\ = -g (n % + ^ x(q) dq) + ^ , dq) - 



x(q) 



where the short-hand notation q ± designates the limit from below/above. Further 
performing integration by parts, and using x(q > qu) = 1, x(q < go) = n, we finally arrive 
at 
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(115) 



f<ik / fid \ 
+ g(Qd-qk) + / g / x(q)dq) 

Jqo \Jq J 



We are actually interested in the replica limit n — > 0. According to the Parisi prescription 
in such a limit the inequality Eq. (11041) should be reversed: 



n = < mi < "^2 < • • • < T^fc < Wfe+i = 1 



;ii6) 



and the function x(q) is now transformed to a non-decreasing function of the variable q 
in the interval qo < q < qk, and satisfying outside that interval the following properties 



x (q < Qo) — 0, and x(q > q^) = 1. 



(117) 



In general, such a function also depends on the increasing sequence of k parameters mi 
described in Eq. (T4"3"l) . 

The form of Eq. (II 15j) makes it easy to perform the limit n — > explicitly, and to 
obtain after exploitation of Eq. (j42|) an important identity Eq. (j44|) helping to evaluate the 
traces in the replica limit. Finally, let us mention the existence of an efficient method 
of the "replica Fourier transform" allowing one to diagonalise (and otherwise work) with 
much more general types of hierarchical matrices, see [23, EH] for more detail. 



B Appendix: Stability analysis of the replica solu- 
tions. 

B.l General relations. 

In this part we are going to derive the general stability equations in the \x— dominated 
regime, and then indicate modifications arising in the R— dominated case. 
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In general, the stability analysis amounts to expanding the function $(Q) in Eq.(14) 
around the extremum point up to the second order in deviations: $ = $sp + 5$ + \5 2 ^. 
In the \i— dominated regime the stationary point is inside the integration region and the 
stationarity condition amounts to 5$ = yielding the system (16)-(17). The term <5 2 <3> is a 
quadratic form in independent fluctuation variables 5q a b, a < b and can be generally writ- 
ten as 5 2 $ = J2(ab),(cd) 5q{ab)G{ab),{cd) 8q(cd)- As usual the stable extremum corresponds to 
the positive definite quadratic form, and along the critical line the quadratic form becomes 
semi-definite. Checking positive definiteness of 5 2 Q amounts to finding the (generalized) 
eigenvalues A of the matrix G^^cd) = dq b 8q d ' s °l v i n g the equations 

X G^ a b) t ( c d)f](cd) = A X C(ab),(cd)V(cd), a < b, (H8) 
(cd),c<d (cd),c<d 

Here the notation 7]^, a < b is used for n(n + l)/2 components of a (generalized) eigen- 
vector rj and C > can be any real symmetric positive definite matrix used to define a 
suitable scalar product (Sq^ A \5q^ B ^ in the space of fluctuation vectors 5q, that is 

(SqW,6qW) = £ bq\tb) G{ab) , (cd) $q{cJ) • (119) 

(ab),(cd) 

Indeed, for any choice of C all generalized eigenvalues are real due to the symmetry prop- 
erties of the matrix G , and the eigenvectors 77^ and 7/^ corresponding to different eigen- 
values Aj 7^ Aj are orthogonal with respect to the chosen scalar product: \r}^ l \rj^J = 0, 
hence linearly independent and form a basis. Expanding the fluctuations in this basis 
as 5q = YUPiV® one then finds generically <5 2 $ = J2j h-iPl c i an d the positive 

definiteness amounts to the condition Aj > for all i. The instability occurs when some of 
the eigenvalues vanish: Aj = 0, and it is easy to see that the condition for the instability is 
independent of the choice of the matrix C in the definition of the scalar product Eq. (I119j) . 
In fact, in our analysis we find it convenient to introduce a formal convention r}( ab ) — V(ba) 
for eigenvector components, which makes it natural to think of the eigenvectors 77 as being 
real symmetric matrices. Accordingly, we find it convenient to define the scalar product 
in that vector space for any two such eigenvectors rf A ' and as 

(yW,^) =Tr[^ A V fl) ] (120) 

This simply corresponds to choosing the matrix C in Eq. (11191) to be diagonal, with the 
diagonal entries given by C( afe ) >(ab ) = (2 - 8 ab ). 

Our first task is to determine explicitly the entries of the stability matrix G for our 
problem. The structure of the replica free energy functional Eq. (JTJJ suggests to represent 
the matrix G as a sum of two terms G = Gi + Gjj such that: 



t d 2 , , . „ , 1 d 2 _ 

in 1 

(a^b) 



(G 'W->-2«™< < G "W-> = "2T« £ < 121 > 
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where the last summation is assumed to go over all possible pairs (ab) with a ^ b, and we 
introduced a short-hand notation D ab = \(q aa + qbb — 2g a ft)- A straightforward application 
of the Wick theorem shows that: 



(Gi 



>(a?b),(c^d) 



T 



(o-'L («-')« +(«-'L(«- 1 ) J. (° 



' (aa),(cc) 



T 
~2 



Q- 



ad 



122) 



In view of the convention rf( ab ) = rfc ba \ for eigenvector components, we have 



T 

E (Gi) iabUcd) vied) = 7-7; E (Q~ l ) ac M h k^. 

i),c<d ^ ° afc (cd) 



(123) 



(124) 



(cd),c<d 



where 5 ab stands for the Kronecker symbol, and the summation in the right-hand side 
goes over all pairs of indices without restrictions. 

Further simple differentiation gives for the entries of Gjj with a 7^ b the expressions: 



1 



1 



(^-f)(o^6),(o^6) ~~ j, f (Dab), (Gll)(pe£b),{aa) ~ (Gil) ( a ^b),(bb) ~ r^j, 



— f"(D ab ), (125) 



and (Gn)(~£b\ ( c d) = f° r a H other choices of (cd). If however a = b, we similarly find 



(Gll)( aa ),(aa) ~ ^rp 



^rp E f"(D ac ), (Gu)( aa ^( ah ) — — f"(D ab ), (Gn)( aa ),(bb) — ^rpf"(D ab ) 



c, c^d 



and (G 



(aa),(cd) 



for all other choices of (cd). We therefore see that 



E ( G n)(ab),(cd) V(cd) = 7^f"(D ab ) \r] (aa) + r] (bb) - 2r] iab) ] , a ^ 6, 



(cd),c<d 

and for a = b 



(126) 



(127) 



E ( G Il)(aa),(cd) V(cd) 
(cd),c<d 



AT 



E f"( D ac) Waa) + V(cc) ~ ^V(.ac) 



c, c^a 



128) 



Combining all these expressions and definitions we see that the system of equations 
Eq. ( 1118}) for eigenvector component rf( ab ) and eigenvalues A of the stability matrix G 
in the fx— dominated regime takes the form: 

T E (^L (^L^M + ^f"( D oh) Waa) + 77(66) - 2 V{ab) } = 2A V(ab) , a b (129) 



(cd) 



2T 



34 



and for the eigenvector components with a = b: 



(cd) 



Q 



-i 



ad 



V(cd) f"( D ac) [V(aa) + V(cc) ~ %V(ac) 



AV(aa)- (130) 



Assuming T > and introducing the notation A* = 2A/T, we can rewrite the above pair 
in the unified form: 



E(«- 1 )„ c "m( < 3- 1 

cd 



db + ^f"(D ab )(5D ab ) 



1 

2^2 



5 ab J2f"(D ac )(5D ac )=A\ ab) (131) 



where we used the notation 



6D 



1 



ab 



V(aa) + V(bb) ~ 2r/( a6) 



(132) 



Let us now indicate minor changes required in the similar procedure for R— dominated 
regime. First of all, in that regime the diagonal entries q aa are fixed to the boundary value 
q aa = R 2 ,\/a. The stability for the linear deviations 5$ > implies d<&/dq aa < at fixed 
q ab on the boundary, which amounts to the inequality 



b 



ab ■ 



133) 



On the other hand, in the \x— dominated regime the last inequality is replaced by equality. 
Hence, it is easy to see that the inequality (11331) means nothing but \i < fio with T = 
Tb(fio). 

As to the quadratic form 5 2 $, the entries of the matrix G involve in this regime only 
pairs of indices with (a ^ b) and (c^ d). The explicit expressions for those entries are 
given by formally the same equations as in the \x— dominated regime. The eigenvectors i] 
must now have only components t]( a <b) ■ It is again convenient to introduce the convention 
V(ab) = V(ba),a 7^ b for those components, together with r]( aa ) = 0,Va. The corresponding 
eigen-equation then amounts to 



(cd) 



Q- 1 



i,d 



V(cd) - Tj^f" ( R2 ~ Qab) V(ab) = A* 77 (a6 ), a ^ b. 



(134) 



Now we proceed to the analysis of the main equations Eq. (11 2911 1301) and (11341) in various 
cases. 



B.2 Stability of the replica symmetric solution. 

The entries of the matrix Q in this case are given by q ab = qo + (qd — qo)8 ab , so its inverse 
Q^ 1 has the same form (Q -1 ) a b = Po + (Pd~ Po)8 a b, with p and pd defined in Eqs. ()19ll20p . 
Taking this fact reduces the equation Eqs. ()129p for a ^ b to the form (see Eq. ()132p ): 



Po J2V(cd) +Po(Pd 

(cd) 



Po) 



V(ad) + E V(bd) 
. d d 



+ (Pd ~ Po) 2 V(ab) + 7j^f"{D)5D ab 



A* r] {a b) 
(135) 
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where D = q d — q . The equation Eg. (11301) for a = b is similarly reduced to 



PO £ V(at) + 2Po(Pd - Po) £ V(ad) + (Pd ~ Po) 2 V(aa) ~ 7f^f"( D ) £ ^Aic = A * Va 
(cd) d 1 c 

(136) 

Now we can follow faithfully the lines of the classical work by De Almeida and Thouless 
[36] and to provide an explicit construction of the eigenvectors with components 77( a &)- 

The first ("longitudinal") family consists of those eigenvectors which have the same 
form as the replica-symmetric Q matrices themselves: tj^L = (3 + (a — /3)5 a b. Then 

Y,dV(ad) = & + (3 {n — I), and Eqs. (ll35fl36l) is reduced in the replica limit n — > to a 
system of two equations 

2po(Pd - Po)(a -/3) + (p d - po) 2 (3 + ^f"(D)(a - (3) = A* (3 (137) 

2po(Pd - Po)(a -13) + ( Pd - Po) 2 a + ^f"(D)(a - (3) = A* a (138) 

One can find the corresponding two eigenvalues easily, but for our goals it is enough 
to notice that for any values of the system parameters the eigenvalue A* can not vanish 
as A* = immediately implies (3 = a = 0. Hence, fluctuations in the longitudinal family 
cannot induce instability of the replica-symmetric solution. 

Second family r/^ of the eigenvectors can be constructed as those characterized by 
one particular replica index singled out among all n. As all these choices are equivalent 
up to permutation of replica indices, we can take for definiteness a — 1. At this special 
choice of the replica index we have explicitly 

^(li) = w > = V(tl) = t, V6 ^ 1, and ^ = /x + 5 ab {u -/x),Va^l;V6^1, (139) 

To ensure orthogonality of such an eigenvector to those of the longitudinal family, i.e. 
(r]( II \r)( I ' > ) = 0, see Eq. (1120j) . amounts to imposing the conditions 

n — 2 

w + (n-l)i/ = 0, r + /i^— = 0, (140) 

It is easy to check that these conditions in fact ensure that 

£i#3 = -(n-l)(i/-r), EffiW-^Vc^l, so that £^ = 0. (141) 

d d (cd) 

Note that in the replica limit n — >• the relations Eq. (11401) imply u = v and r — \i. 
Then it is easy to check that the equations Eq. (I135II136P for a = 1 are reduced to precisely 
the same system of two equations as that in Eq. fll37lll38|) . with the correspondence a — > 
v, f3 — ► t. Moreover, for any a > 1 the equations Eq. (11351)1361) yield once more the 
same pair of equations. We therefore conclude that the fluctuations corresponding to this 
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second family of eigenvectors are also unable to induce instability of the replica-symmetric 
solution. 

Finally, the third family of eigenvectors r]^ ni > orthogonal to the first two turns out to 
be characterised by a chosen pair of nonequal replica indices (a, b), a ^ b. For the sake 
of defmiteness we can take a = 1, b = 2 when the corresponding eigenvector components 
are given explicitly by: 



£1? = 0, Va; n^iP 



V(aa) 



J, Va; nfa' =^ ( 2i / ) /) 



'/(lb) — 'l(2b) — 7(61) — 7(62) 



V>, V6 > 2, 



</) = (l-^)P ; Va>2;V6>2 
where the parameters £,ip,p satisfy the constraints: 

f + (n- 2)^ = 0, 2^ + p(n-3)=0, 

It is easy to check that these conditions in fact ensure that 

£i 7 ? = 0>Va 



(142) 



(143) 



(144) 



so that automatically J2( c d) V( c d) = ®- We t nen nn d that the equations Eq. 
satisfied identically, whereas the equations Eq. fll35j) are reduced to the form 



are 



n (III) 

'l(ab) 



(145) 



As the replica symmetric solution in the //—dominated regime implies D = T/fi, Pd~Po = 
p/T, see Eqs. (l2~Tl |24"|) . we infer from Eq. (I145j) that the stability condition A* > indeed 
implies the inequality Eq. flBTT) . 

Turning to the stability of R— dominated replica-symmetric solution, we find that the 
eigen-equations Eq. (I134j) are reduced to 



Po ^2v(cd)+Po(Pd-Po) 

(cd) 



V(ad) + V(bd) 



L d 



+ (Pd-Po) 2 V(a 



(146) 



-Tjvf" ( r2 ~ ?o) V(ab) = A* T] {ab) , a^b 



J^2 

The subsequent analysis is completely analogous to one performed above in the //—dominated 
regime. In particular, the relevant eigenvalue again corresponds to the third family of 
eigenvectors satisfying the constraint Eq. (jl44p . This fact, together with the relation 
Pd — Po = l/(-R 2 — Qo) immediately yields the inequality Eq. f[3"2l as the corresponding 
stability condition. 
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B.3 Stability of the Parisi solution with one step of the replica 
symmetry breaking. 

In the following we will show that in this situation there are essentially two relevant 
eigenvalues, both in the \x— dominated and R— dominated regimes: 

= 7 ~ 7^f"(qd - Qi) (147) 

[Qd - qi) T 2 

and 

K = — ^ ^ - ^f"(q d ~ qo) ■ (148) 

If both are positive, all other eigenvalues are positive and the system is stable. In the 
replica symmetric case both eigenvalues fall together. Let us in general note that if we 
know the solution on the line T = !&(//), which is defined by the condition q d = R 2 , we 
know the solution in the whole i?-dominated regime, since there the values qj , q\ , qo and m 
are independent of fi (see the main text of the paper). Thus it is in general enough to 
proof the stability in the /^-dominated regime. 

To this end, let us consider the function R(q) = S(q)/T, where S(q) was defined in 
Eq. (]4"9]) of the main text. Specified for the case of 1RSB this function satisfies (cf.(|69l)): 



R(q ) = 0, R{q x ) = 0, and R(q) dq = . (149) 



Moreover we know that R(q) must change sign once between q and q±. Let us calculate 
the positions of the corresponding two extrema g^ 1 ' 2 -* from R'(q e ) = 0: 

1 1 f"{qd ~q e ) = 0. (150) 



(qd - qi + m(q 1 - q e )) 2 T 2 
which yields 



T 

q d ~qi+ m(q 1 - q e ) = = (151) 



for q e = g^ 1 ' 2 -*. Since the right-hand side is convex according to our assumption for short 
range correlations, we have g^ — gi + m(qi — g) > T / J f"(qd — g) for q^ < q < q( 2 \ 
hence R'(q) < in that interval. This however means that for g ^ (?e j 9e ) necesserily 
R'(q) > 0, in particular R'(qo) > and R'(qx) > 0. (The accidental case R'(qo) = or 
R'(qx) = happens on the AT-line). Note however that R'(qo) is exactly equal to Aq and 
similarly R'(qi) = A* K . This observation proves stability everywhere in the \x— dominated 
phase. The stability condition becomes marginal along the AT-line. 

In the rest of this appendix we will sketch the derivation of the fluctuation eigenvalues 
in both regimes, /z-dominated and -R-dominated. First let us define a Parisi-block matrix 
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P ab = 1 inside a diagonal Parisi m x m block and zero outside. With this definition the 
entries of the Parisi matrix Q has the form 

Qab = Qo + (qi ~ qo)Pab + (qd ~ 9l)*crf> • ( 152 ) 

The inverse matrix Q _1 has the same form: 

(Q _1 ) aft = A + (A, - A )P ab + (A d - A 1 )S ab . (153) 

The matrix P ab can be used to introduce the following shorthand notations for several 
types of averages of the eigenvector components r\i a b) over the indices inside a Parisi block: 

11 1 

Vab = — (V P )ab, Vab = ~^{PV P )ab, Vaa = — ^2V(cc) Pea • (154) 

With these definitions the eigen-equations Eq. (11311) can be written as 
A*7] ab = A 2 J2 V(cd) + {A x - A ) 2 m 2 7]- al) + (A d - A^rj^ 

c.d 

+mA (A l - A ) J2(riac + Vdb) + M^d ~ A x ) J2(V(ac) + r/ (c6) ) (155) 

c c 

+(Ax - A ){A d - Aa)m(77a 6 + rj al ) 
+^f"(D ab )5D ab - 5 ab £ ^f'\D ac )5D ac . 

Here, explicitly 

f"(D ab ) = f"(q d - g ) + (f"(q d - Ql ) - f"(q d - q ))P ab (156) 

and 

2 f"(Dac)SD ac = f"(q d - g ) In r] (aa) + ^ r/ (cc) - 2 ^ V(ac) I + 

c \ c c / 

+ m [f"(q d - gi) - f"(q d - %)] (v(aa) +Vaa- ^Vae) • (157) 



Now one observes that one can derive a closed system of three equations for the 
following three combinations 

X/*7M> ^2 Vaa, 22V(ab) ■ 
a a ab 

This yields three eigenvalues in the limit n — > 0: 
A\ = (A d -A 1 + m(A 1 -A )) 2 

I 1 TTi 

^ = (A d -A i y-—f"{q d - qi ) + ^ Fr {f'{ qd -q 1 )-f{ qd -q Q )) 

A; = (A d -A 1 +m(A 1 -A )f-^f"(q d -q Q ) . (158) 
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As the next family of eigenvectors we use those satisfying the constraint J^aViaa) — 
J2aVaa = J^abViab) = 0. One can show that constraints of this sort ensure orthogonality 
of the new families of eigenvectors to the old one (compare with the de-Almeida-Thouless 
analysis of the previous section). Subsequently, one derives a system of equations for 

b 

In the limit n — > that system yields the same eigenvalues A{, A%, A3. 

Using now as the constraint the conditions rfaa = Vaa = J2b Vab = simultaneously 
for all a we derive a single equation for the component with the indices a and b in 
different Parisi blocks (we write a^b). The resulting eigenvalue is again A3. 

In the next step we impose another constraint %^ = for all a, b to obtain a set of 
equations for 

V(aa), Vaa, ^2V(ab) ■ 
b 

This leads to another three eigenvalues 

A* = (A d - A^Aa- A 1 + m(A 1 - A Q )) 

1 9 _ m 

A5 = (A d -A 1 y--f( qd - qi ) + ^ r (f( qd - qi )-f( qd - qo )) 

A* 6 = (A d -A 1 )(A d -A 1 + m(A 1 -A ))-^f( qd - qo ). (159) 

We then obtain a single equation for 77^, which leads for a 7^ b to Ag again. The 
remaining equation for rj^ yields for a 7^ b a new eigenvalue 

A* 7 = (A d -A 1 ) 2 -^f"( qd - % ) (160) 



and for a = b and a 7^ b 



A* 8 = (A d - Ax) 2 - ±f"(q d - qi ) . (161) 



Assuming f"(x) monotonously decreasing, < m < 1 and < go < Qi < Qd the relevant 
eigenvalues are obviously A^ = Ag and Ag = A^. Here we use that the eigenvalues of the 
matrix Q are 

Qd-Qi, qd-qi + m( qi -qo), q d - qi + m( qi - q ) + nq (162) 
These are inverses of the eigenvalues of the matrix Q^ 1 : 

A d -A u A d - A 1 + m(A 1 - A Q ), A d -A 1 + m(A 1 - A ) + nA . (163) 
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Turning now to the i?-dominated regime, the eigenvalue equations (11341) with the use 
of the short-hand notations introduced in Eq. (11321) can be written in the form: 



A*5D ab = (A 2 Q ^2^D c , d + (A 1 -A ) 2 m 2 5D- a - b +(A d -A 1 ) 2 5D ab 

c,d 

+ mA (A 1 - Aq) ^D ac + 5Q cb ) + A (A d - A 1 ) + 5D cb ) 



+ (At - A )(A d - A x )m{8D- ab + 5D ab )){\ - 6 ab ) - ^f"(D ab )5D ab 
We proceed in the same way as before. The equations for 



produce the quadratic equation 



a ab 

A- A* B 
C D-A* 



with 



.4 
B 



(m - l)(m(A 1 - A ) 2 + 2{A X - A )(A d - Ax)) + {A d - Axf 



J^2 



f"{q d - Qi) 



(m - l){nA 2 + 2A (A d -A x + m{Ax - Aq))) 



C = - m (A 1 -A ) 2 -2(A 1 -A )(A d -A 1 ] 



J^2 



f(q d - Ql ) + ^f"(q d - g ) 



D = n(n-l)A 2 + 2(n-l)A (A d -A 1 +m(A 1 -A )) + (A d -A 1 +m(A 1 -A )f 
-TpfifLd ~ ?o) ■ 



In the next step we obtain for 



5D- n 



b 



another quadratic equation 



A - A* B 
C D-A* 



with 



1 



A = ( m -l)( m (A 1 -A y + 2(A 1 -A )(A d -A 1 )) + (A d -A i y-—f"(q d -q 1 
B = ( m -l)(2A (A d -Ax + m(Ax-A ))) 



J"2 ' 



C = - m ( Al -A ) 2 -2(A 1 -A )(A d -A 1 )-—f"(q d -q 1 ) + —f"(q d -q ) 



J>2 ' 



J^2 • 



D 



(n - 2)A (A d — Ax + m(Ax - A )) + (A d -A x + m(Ax - A )) 2 - j^f"(q d - q Q ) 

Both equations fall together for n — > and produce the same eigenvalues A^, A^. One 
can show under the conditions we have that both eigenvalues are positive provided Aq and 
A^ are. Next one derives a single equation for SD^ b for a ^ b leading to the eigenvalue 



A* = (A d -Ax + m(Ax - A,)) 2 - ±f(q d - q ) 



(164) 
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Note that A3 here is formally equal to A3 found in the previous analysis of the /1— dominated 
regime. On the next lower level we obtain for 

b 



the eigenvalue equation 



A - A* B 
C D-A* 



with 



A = ( m -2){A l -A )(A d -A 1 )) + (A d - A^ 2 - ^f"{q d - qi) 

B = (m — 2)A (A d — A K ) 

C = - 2 (A 1 -A )(A d -A 1 )-^f"(q d -q l ) + ^f"(q d -q ) 

D = (A d -A 1 )((n-2)A + m(A 1 -A )) + (A d -A 1 ) 2 -^f(q d -q ) 

leading to A4, A£. Again both are positive if Aq and A^ are. On the next lower levels 
we again reproduce the previous eigenvalues Ag, Ay and Ag of the \i— dominated analysis. 
As the result in both regimes are Aq and A* K the relevant eigenvalues, which have to be 
positive for stability - and they are indeed positive as we have demonstrated before. 
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Figure 2: The phase diagrams for the potential with logarithmic correlations, with g = 
2, a = y/2. In the case (a) the radius of the box is R — 1 < R cr = y/2, in the case 
(b) R = \fh > R cr , in the case (c) R = oo. The case (d) corresponds to the choice of 
parameter \x = 0. The correct boundary Tj, between the R— dominated and \i— dominated 
glassy phases for the case (b) is represented by full vertical line. The notation for phases 
are the same as in fig.l. 
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Figure 3: The phase diagrams for the potential with short-range correlations, f(x) = e~ x . 
In the ) the radius of the box is R = 1 = R cr , in the case (b) R = V3 = R*. Dotted 

line represents the wrong branch of the boundary T b between the R— dominated and //— 
dominated glassy phases and should be replaced by the full line close to vertical. The 
notation for phases are the same as in fig. 1. 
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Figure 4: The phase diagrams for the potential with short-range correlations, f(x) = e~ x . 
In the case (a) the radius of the box is R = > in the case (b) R — oo and the 
case (c) corresponds to /i — 0. The dotted curves in (a)-(c) represents the transition line 
Ti(/i) found from the condition m — > 1 and the broken curve in (B) and (c) represents 
the part of Tat line below Ti. The notation for phases are the same as in fig. 1. 
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